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The HERAFitter project provides a framework for the determination of parton distribution 
functions (PDFs), and tools for assessing the impact of new data on PDFs. In this contri¬ 
bution, HERAFitter is used for a QCD analysis of the legacy measurements of the W-boson 
charge asymmetry and of the Z-boson production cross sections, performed at the Tevatron 
collider in Run II by the DO and CDF collaborations. The Tevatron measurements are included 
in a PDF fit performed at next-to-leading order, and compared to the predictions obtained 
using other PDF sets from different groups. The measurements are in good agreement with 
NLO QCD theoretical predictions. The Tevatron data provide significant constraints on the 
d-valence quark distribution. 


1 Introduction 

According to the factorisation theorem, cross sections in hadron collisions are calculated by con- 
voluting short distance partonic reactions with parton distribution functions (PDFs). Discovery 
of physics beyond the Standard Model at hadron colliders relies on the precise knowledge of 
the proton structure. Moreover, PDFs are among the dominant uncertainties for the measure¬ 
ment of the W mass, and for gg —>• H production. HERAFitter 1 provides a framework for the 
investigation of various methodologies in PDF fits, and tools for assessing the impact of new 
data on PDFs. It is widely used by the LHC experiments to improve the sensitivity of new 
measurements to PDFs. Full information about the project, downloads and documentation can 
be found at herafitter.org. A schematic view of a PDF fit, as implemented in HERAFitter, is 
shown in Fig. [l] 

In this contribution, HERAFitter is used to perform a QCD analysis of the Tevatron Run II 
legacy measurements of the VF-boson charge asymmetry and of the Z-boson production cross 
sections 2 . At the Tevatron proton-antiproton collider, the production of W and Z bosons is 
dominated by valence-quark interactions. Whereas the primary source of information on the 
proton PDFs comes from deep-inelastic scattering (DIS), Drell-Yan production of W and Z 
bosons in proton-antiproton collisions can provide additional information, particularly on the 
d-valence quark PDFs. 
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Figure 1: Schematic representation of a PDF fit in HERAFitter. 

2 Tevatron W and Z measurements and QCD settings 

The most recent measurements of VF-boson charge asymmetry and Z-boson inclusive produc¬ 
tion performed in Run II of the Tevatron collider are considered. They include the Z-boson 
differential cross section as a function of rapidity, measured by DO 3 } the Z-boson differential 
cross section as a function of rapidity, measured by CDF ® the charge asymmetry of muons as 
a function of rapidity in W —> fiv decays, measured by DO 5 } the VF-boson charge asymmetry 
as a function of rapidity, measured by CDF 6 } the VF-boson charge asymmetry as a function of 
rapidity, measured by DO '. Besides the Tevatron VF- and Z-boson measurements, the HERA 
I combined measurements of the inclusive DIS neutral- and charged-current cross sections mea¬ 
sured by the HI and ZEUS experiments® are used. 

In general, the correlation model of the experimental uncertainties recommended by the 
Tevatron experiments is adapted and followed in the QCD analysis, with the exception of the 
experimental systematic uncertainties related to trigger and lepton identification efficiencies, 
which are treated as uncorrelated bin-to-bin. 

The QCD analysis and PDF extraction is performed with the open-source framework HERA¬ 
Fitter. The charm mass is set to m c = 1.38 GeV, as estimated from HERA charm production 
cross section®, and the bottom mass to mb = 4.75 GeV. The strong-interaction coupling constant 
at the Z boson mass, a s (Mz), is set to 0.118, and two-loop order is used for the running of a s . 

The PDFs for the gluon, u-valence, d-valence, u, d quark densities are parametrised at the 
input scale of Qq = 1.7 GeV 2 . The contribution of the s-quark density is taken to be proportional 
to the d-quark density by setting xs(x) = r s xd(x), with r s = 1.0. The strange and anti-strange 
quark densities are taken to be equal: xs(x) = xs(x). 

The impact of a new data set on a given PDF set can be quantitatively estimated with a 
profiling procedure 10 . The profiling is performed using a \' 2 function which includes both the 
experimental uncertainties and the theoretical uncertainties arising from PDF variations: 
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The correlated experimental and theoretical uncertainties are included using the nuisance pa¬ 
rameter vectors /3 ex p and /3th, respectively. Their influence on the data and theory predictions 
is described by the TU 03 and T 16 matrices. The index i runs over all ZVdata data points, whereas 
the index j (k) corresponds to the experimental (theoretical) uncertainty nuisance parameters. 























Table 1: Results of a 15-parameters fit to the to the HERA I and Tevatron W- and Z-boson 
data. The contribution to the total Xmin °f each data set and the corresponding number of 
points are shown. 


Data set 

X 2 / number of points 

HERA I 

515 / 550 

DO da(Z)/dy 

23 / 28 

CDF da{Z)/dy 

33 / 28 

DO muon charge asymmetry in W —> 

12 / 10 

CDF W charge asymmetry in W —> ev 

15 / 13 

DO W charge asymmetry in W — > ev 

16 / 14 

Total Xmin / dof 

615 / 628 


The measurements and the uncorrelated experimental uncertainties are given by <r ? cxp and Aj, 
respectively, and the theory predictions are aj h . 

3 Results 

A QCD fit analysis performed on the Tevatron W- and Z-boson data, together with the HERA 
I data, is used to assess the impact of the Tevatron data on PDFs. The profiling is used to 
assess the impact of the Tevatron data on various PDF sets. 

The optimal functional form for the PDF fit, which corresponds to 15 free parameters, is 
found through a parametrisation scan, and is used for a fit to the HERA I data only, and for a 
fit to the HERA I and Tevatron W- and Z-boson data. Table [l] shows the Xmin P er degrees of 
freedom (dof) of the fit to the HERA I and Tevatron W- and Z-boson data. The contribution 
to the total Xmin °f eac h data set, referred to as partial x 2 , is also shown. The partial % 2 per 
number of points of each of the Tevatron and HERA I data set is close to unity. 

The central value and the uncertainties of the PDFs are evaluated with MC replicas®. 
Figure [2] shows the comparison of the PDFs extracted with the MC replica method by fitting 
the HERA I data, and by fitting the HERA I and Tevatron W- and Z-boson data. A significant 
reduction of the PDF uncertainties is observed in the fit which includes the Tevatron W- and 
Z-boson measurements, in particular for the d-valence quark. 

The impact of the Tevatron W- and Z-boson measurements on the CTlOnlo 12 and MMHT2014 13 
is assessed by profiling. The uncertainties of the CTlOnlo PDFs are scaled to 68% confidence 
limit. The compatibility of the Tevatron data with the CTlOnlo, MMHT2014 and NNPDF3.0^I 
sets is tested by evaluating the x 2 function of Eq. ([!]). The partial x 2 per number of points of 
each of the Tevatron data set, and the total x 2 / dof, are close to unity for all the PDFs. 

The CTlOnlo and MMHT2014 PDFs are profiled to the Tevatron W- and Z-boson data. 
The results of the profiling on the relative uncertainty of the d -valence PDF are shown in Fig. [3| 
Significant reduction of the uncertainty is observed for both sets. 

Tables of the Tevatron measurements, with updated correlation model, and corresponding 
APPLGRID theoretical predictions 15 are publicly available at herafitter.org. 
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Figure 2: (a) PDFs at the starting scale Q 2 

determined with a fit to the HERA I data (blr 
W- and Z-boson data (yellow), (b) Relative PI 



(b) 
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3 ), and with a fit to the HERA I and Tevatron 
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Figure 3: Relative uncertainties of the d-valence PDF at the scale Q 2 = 1.7 GeV 2 as a function 
of Bjorken-x before and after profiling for the (a) CTlOnlo and (b) MMHT2014 PDFs. 
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